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Abstract : 

We elaborate on an observation of Maes- van Wieren (2006) to obtain fluctuation 
symmetries also for time-symmetric quantities. Examples are given, analytic 
and numerical, yielding time-symmetric path-observables with fluctuations sat- 
isfying a Gallavotti-Cohen type symmetry. From these results one is actually 
introduced to stationary nonequilibrium by a different phenomenology. It deals 
with a complementary class of what we may call active fluctuation symmetries, 
again general non-perturbative nonequilibrium relations but not expressed in 
terms of the traditional dissipative variables; they rather involve the notion of 
dynamical activity. In particular, we derive Green-Kubo like relations for differ- 
ences in dynamical activity. The illustrations include boundary driven Kawasaki 
and zero range models and the spinning Lorentz gas. 



1 Introduction 

General and non-perturbative relations are not so common in nonequilibrium 
physics. Recent decades have therefore seen a big interest in the fluctuation sym- 
metries of the entropy production as pioneered in the papers [1] [2] . It was found 
that these symmetries are an expression of local detailed balance, implying that 
the total path-wise entropy flux is the source term of time-reversal breaking in 
the nonequilibrium action; see [3l |4j [5l IH [7] . Local detailed balance refers to the 
underlying microscopic time-reversibility that governs the contact between the 
system and each (equilibrium) reservoir in the environment [3 [HI HI UHl [HI [H] • 
Also the nonequilibrium free energy relations, called Jarzynski relation after 
[13j . are of a very similar nature. 

The natural place to study all these fluctuation relations is in large devia- 
tion theory for occupations and currents in nonequilibrium systems. The present 
paper wants however to take some distance from the original work which con- 
centrated mostly on heat and entropy production. Here we add fluctuation sym- 
metries for aspects of dynamical activity, which belong to the time-symmetric 
fluctuation sector of a nonequilibrium system. This dynamical activity (or frc- 
nesy, as called in the context of linear response |14] ) is an important kinetic 
aspect of nonequilibrium and is of growing importance in fluctuation theory. 
We refer to [15l[T6| for some impressions and further references. The main re- 
sult of the paper (in Section [4|) gives fluctuation symmetries for differences in 
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dynamical activity. 

For the plan of the paper, the next section repeats the central idea after 
[T71 131 HH m] • Besides time-reversal symmetry we add a second symmetry which 
can be spatial or internal and that gives rise to additional fluctuation symme- 
tries. The standard example of a fluctuation symmetry for the entropy flux is 
reviewed in Section |3l We mention also how we view the relation with other 
topics in nonequilibrium statistical mechanics, in particular and less known, 
how the Kubo formula (and not just the Green-Kubo relations) follow from 
the fluctuation symmetry. Staying still in the same context we then treat two 
examples (boundary driven Kawasaki and zero range dynamics) in Section 2] 
that are explicitly calculated to give fluctuation symmetries for differences in 
their dynamical activity. There we find the main point of the paper, giving 
fluctuation-activity relations. Section [5] applies that to the spinning Lorentz 
gas where the notion of dynamical activity gets a specific physics realization. 
Computer simulations validate there our guesses also in the non-Gaussian fiuc- 
tuation sector. 

2 General observation 

Wc start with the formal content of fluctuation symmetries, [T71 [5]. 
Let X € denote a fluctuating quantity. That means that its outcome is 
variable and uncertain as, in physical terms, it depends on hidden or more 
microscopic degrees of freedom. For mathematical modeling we will consider 
probability distributions for X where fl is the space of possible outcomes (in the 
universe). This in turn also means that D, supports some elementary structure 
such as used for integration. More relevant is the presence of involutions Q 
and r on 51, preserving there such elementary structure, and being mutually 
commuting, = 9^ = Id, OF = FO. 

We will always have a reference probability law Pq for X which is both Q— and 
F— invariant; Po{QX) = Po{X) = Po{TX). Our main interest is in a probability 
law P on X which we assume has a density with respect to PqI 

dP(X) e-^(^MPo(X) (2.1) 

for "action" A, which for our purposes is mostly explicitly known. For all the 
meaning of (|2.ip we pretend for a moment that X takes a finite number of 
values so that expectations (•) under P are simply written as finite sums 

{fix)) = J2 f{^)p{x) = E /(^) e"^^'^ Poi^) 

X X 

for an observable / and with P{x) the probability that X ^ x. 
Examples 

1. Dynamical ensembles. $7 can be the set of allowed (possibly coarse- 
grained) trajectories of a dynamical system. The trajectory X is then a sequence 
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of states {xt,t £ [0, T]) of a (sub)systein over a time-interval [0, T], and we could 
take = time-reversal, i.e., &X = {TTXT-t,t G [0,T]) with tt the kinematical 
time-reversal such as flipping the momenta of all particles. 
As an illustration of a possible origin of dynamical ensembles we take a collection 
of sites i G V where we think of V as the volume and to which we assign coupled 
oscillators {qi,Pi) £ M^. The dynamics is Hamiltonian with potential U except 
at the boundary sites i G dV where we add Langevin forces: 

dq, = p, dt, leV (2.2) 
dU 

dp, = - — {q)dt, ieV\dV 
aqi 

dp., = -^iq)dt--/n,p,dt+ J'^dW.it), i£dV 
oqi V Pi 

The Wi{t) are independent standard Brownian motion. When Kj/J^ ~ /? for all 
i € dV, then p*^ cx exp —f3{p^ /2 + U{q)} is a stationary equilibrium distribu- 
tion for all 7 > 0. The stationary process for (|2.2p is then reversible, giving 
the reference Po for which Po = PoQ- For trajectory X = {{pt,qt),t G [0,r]) 
we have the time-reversed trajectory {QX)t — {—pT-t,qT-t)- When = 1 
in ()2.2|) appears the nonequilibrium process P which can either be started at 
time zero from or from its own stationary density making then a stationary 
nonequilibrium process. When (always for Ki = 1) the inverse temperatures 
Pi of the reservoirs connected to dV are different, we break time-reversal in- 
variance and also spatially some symmetry will be broken (e.g. between left 
and right when ^ is a linear chain) even when the potential U is homogeneous. 
The second symmetry F can then for example reflect the volume (reversing left 
and right, which is FX = {pt{L — i),qt{L ~ i),t G [0,T],i = l,2,...,i) for 
V = {0,1,2..., L},dV = {0, L}. Note of course that P{TX) =^ P{eX). 
The action A for (|2.1[) was calculated in Section 3.1 of pUj . 



2. Random matrices Consider the non-centered Wishart random matrix 
model; see e.g [21] for the statistic and signal processing literature. Here X is 
axL N X N complex matrix with distribution 

dP(X) = 1 e-^Tr[(X-B)-(X-B)] 

where dX = Y[dXijdXij stands for the volume element in . The matrix 
B is given, and its presence breaks the unitary invariance QX = U X U* for a 
unitary matrix U. The reference ensemble has B = 0, with distribution Pq, so 
that for 

dP{X) = dP„(X)e^Tr[i3*M-s]+c 
for some normalization C that depends on B and N. 

3. Binomial distribution Let X be the number of successes of a repeated 
and independent Bernoulli experiment with parameter a G [0,1]; i.e., X = 
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vi + V2 + +Vn with Vi = 0, 1 independently with Prob[wi = 1] = a. Take 
involution QX = n — X, which leaves invariant the distribution Po on X for 
a = 1/2. Then, for a; = 0, 1, . . . , n, 

P{x) = P„(a;)2V(l-p)"--=P„(X)(2(l-p))"(^)" 

P(x) = Pin - x) (_P_f^-'' ^ p(ex) cxp{(2x ~ n) log ^-} (2.3) 
1 — p 1 — P 

is the distribution of X for parameter value a = p. 

We continue with our general observation, starting from (|2.ip . Define 

5 := AQ~A 

T := Ae + A (2.4) 
R yier - A 

so that 

Note also that 2R = 7T — T + S + ST so that i? is the antisymmetric part of 
T under F when 5* is antisymmetric under F; 

sr = -S ^ R = ^{Tr -T) (2.5) 
The very definitions (|2.4p imply the identities 

(/(eX)) = Po(x) = (/(X)e-^(^)) (2.6) 

ifieVX)) = ^ /(x)e-^(«i^-) F„(x) = (/(X)e-«(^)) (2.7) 

X 

for all functions / on O. From (|2.7p we also have for 0— symmetric observables 
/ = /e such as / = TF - r that 

(/r) = (/e-^(^^-'^)-^(^+^^)) (2.8) 

which we call an active fluctuation symmetry for reasons that will become clear 
in Section m 

There is a rewriting to the more familiar Gallavotti-Cohcn type fluctuation 
symmetries. From (|2.6p by taking f{x) = S{S{x) — a) we get 

Prob[S'(X) = -a] = e-" Prob[S'(X) = a] (2.9) 

and, from (|2.7p by choosing f{x) = S{R{x) — r), 

Prob[i?(X) = -r] = e-"" Prob[i?(X) = r] (2.10) 
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with probabilities referring to the probability law P. 

These relations are general and can be applied in a variety of ways. For example, 
as will be seen in Section and in equations p.l4p - p.l6p . from (|2.6p follows 
that for all functions g on X, 

to first order in the action A of (|2.ip and where the last expectation (•)° is with 
respect to Po- That is a (generalized) Kubo or fluctuation-dissipation relation, 
valid around equilibrium Pq. Higher orders and nonlinear response around equi- 
librium can be obtained from combining (|2.6p and (|2.7|) as done in But 
the same holds replacing Q — >• QT and S* — > i? in which case wc will arrive at 
fluctuation-activity relations |23j . 

Another easy consequence is that always {R{X)) > 0, (TT — T) > and 
{S{X)) > 0, for example useful to determine the direction of currents, [19]. 
All these consequences remain basically intact also for variables that differ from 
iS* or i? by a total (time-)difference as long as some boundedness of these terms 
can be assured (referring to the first example above). We then get asymptotic 
fluctuation symmetries, where (|2.6p - (|2.10p are not exact but only valid in some 
limit (of observation time). 

It goes without saying that the relevance of the fluctuation identities (j2.6p - (|2.7p 
depends crucially on the systematic and operational meaning of S and T. It 
was understood before that S is deeply related to changes in entropy (as we will 
repeat in the next Section) , and Sections |3HS] treat examples where T is made 
visible and related to the dynamical activity. 

3 Standard example: entropy flux 

The present section contains the standard application of (j2.6|) to obtain a fluc- 
tuation symmetry for the total entropy flux in a model of nonequilibrium. There 
will be nothing particularly new here, except for the style of presentation, with 
its emphasis on (|2.6|) and some reflections towards the end of the section connect- 
ing the fluctuation symmetry also with response theory. Examples for spatially 
extended systems are not so common in the literature on fluctuation symmetries 
and Section 2] will provide some. 

To be more specific, consider a Markov jump process on a finite state space 
K. We specify the transition rates kt{x,y) (time-dependent) for jumps x — > y 
as 

kt{x,y) = ip{x,y) exp{^[U{x,at) - U{y,at) + F{x,y)]} (3.1) 

where at is a time-dependent (external) protocol changing the energy function 
U. Wc take the driving F{x,y) = -~F{y,x) antisymmetric and the reactivities 
tp{x,y) = ip{y,x) symmetric. The other time-dependent parameter /3t > is 
the changing inverse temperature of the environment (in units where Ub = 1); 
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the nonequilibrium driving sits entirely in the function F and in the time- 
dependence of both the protocol at and the inverse temperature f3f We call 
U{x,a) the energy of the system when in state x for external value a, because 
we imagine that the changes in U are exactly balanced by the change of energy 
in the environment. Clearly, if = and when at = a, pt = P are constant, then 
the process is reversible with stationary distribution p^{x) oc exp — {x , a) . 

For a given path X = {xt,t ^ [0,T]) over the time- interval [0, T] the energy 
changes 

U{xT,aT) - U{xo,ao) = [U{xs,as) - U{xs- ,as)] + / —{xt,at)atdt 

(3.2) 

because of two effects: for fixed value at the state changes and energy is ex- 
changed with the environment as heat 

Qo{X) := [U{xs,as) - U{x,- , a,)] (3.3) 

s<T 

(sum over jump times in X). Secondly, for fixed state Xt the external value 
changes at = doing work 



Wo{X) := 




Equation (j3.2p mimics the first law of thermodynamics. The change in energy 
of the system equals the change in internal energy received as heat Qo from the 
environment plus the amount of work Wo done on the system by the environ- 
ment: 

U{XT, ar) - U{xo, ao) = Qo{X) + WoiX) 

The nonequilibrium driving F can be added and subtracted from that balance. 
We think of it as doing work on the system, which is instantaneously released 
as heat, so that now U{xT,aT) — U{xo,ao) = Q{X) + W{X), but with 

Q(X) :=Qo(X)-^F(xt-,xO, W{X):^Wo{X) + Y,F{xt^,xt) (3.4) 
t t 

all depending on a specific path X. We refer to [53] for more details and insights 
on stochastic energetics. 

In the same spirit we can also associate a change in entropy of the environment 
to a trajectory X . The idea is that the environment consists of big equilibrium 
reservoirs undergoing only reversible changes in interaction with the system. 
We look back at and ((H^ to define 

5o„,(X) := - ^ PsSQs = ^ A {F{xs- , x,)~[U{xs,a,)~U{x,- , a,)] } (3.5) 

S S 
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for the change of the entropy in the environment (always per fc^). That is called 
the entropy flux, which can be split into a reversible part, due to the energy 
exchange, and an irreversible part 

a{X):^Y,PtF{xt-,xt) (3.6) 
t 

We now repeat the observation of [31 1^ that the entropy flux (|3.5p can be 
obtained as for (|2.6p as the source term of time-reversal breaking. 

For a good moment let us leave out the kinematical time-reversal t: on K and 
proceed with the undecorated time-reversal Q defined on paths X via {&X)t = 
Xr-t for t G [0,T]. We check from ([Xi]) that S^^^{X) = -S^^^{&X), or, 
the entropy flux per path is antisymmetric under time-reversal. Let now be 
the path distribution when we start at time zero from the law /i. The time- 
dependence of the protocol can be reversed to define kt{x,y) :— kT-t{x,y). 
Choosing a law ly on K the latter Markov process has a distribution on the 
paths X = {xt,t e [0,T]). Assuming ^,1/ > and (the dynamical reversibility) 
that kt{x^ y) ^ implies fct(?/, x) = 0, we can find the S in (|2.6p by writing 

"^^^ -e^ (3.7) 



dP^e 

and get 

K^o) , hi{xii,xt^)kt^{xt^,xt^) . . .kt^{xt„_^,XT) 



S{X) - log 



v{xt) kt^{xT,Xt„_J ■ . . fcti(xt2,a;tJfctj(xti,a;o) 



for jump times ti,t2, • ■ • in X. Indeed the jump times in &X are respec- 
tively T — tn, . ■ . ,T ~ t2,T ~ ti, and we have substituted kT-t„ixT,Xt„_i) = 
ktn{xT , Xt^_-y) etc. One can see what this becomes for p.ip . Substituting into 
the previous formula makes 

S{X) - log 4^ = E Pt{U{xt- , at) - U{xt,at) + F{xt^ , xt)} (3.8) 

which is p.Sp . That relation can be called a (generalized) Crooks relation [5], 
and for P = it almost immediately produces Jarzynski identities which are 
used to evaluate equilibrium free energies from the fiuctuations of the dissipative 
work — we refer to the literature and the references therein for more details, 
[Ml [251 [7]. 

Let us now specify to the case where /3( = /3, = a are constant in time. 
In particular, with respect to p.6p . and for state functions h^^{x) := log^(a;) -|- 
l3U{x),hu{x) := log u{x) + (iU{x), we have the identity 

S{X)^ [3Y,P(.^t-,xt) + h^{xo)-K{xT) (3.9) 
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for all trajectories X. Note that the left-hand side is defined from p.7p imple- 
menting (|2.6p . while the right-hand side is defined from the heat and (|3.5p - p.6|) . 
Therefore, the identities (|3.8p - p.9p are the core of what is generally called the 
fluctuation symmetry, the fluctuation relations or the fluctuation theorem (tran- 
sient or steady state) for the entropy production [7]. 

3.1 Exact fluctuation symmetry 

In the following we restrict ourselves to time-homogeneous Markov processes 
and we do no longer write the dependence on at = a. We take also /3 = 1. 

Consider the reference reversible process Pq started in equilibrium po for 
which there is detailed balance with rates 

fco(x,y) = V(^,y)e^[^(^)-^(^)1, Po{x) = le-^(^). 

The nonequilibrium process has rates k(x,y) = ko(x,y) expi^(a;, y)/2 and we 
choose to start it also from po- Its distribution on paths X in the time- interval 
[0, T] is then denoted by P. We proceed as in (|2.4p to find 

S{X)=y^F{xt-,xt), T{X) = 2[ [^{x,) ~ Ux,)]ds (3.10) 
t ■J" 

for escape rates ^{x) := J2y ^{^^ v)- ^O"^ clearly ()2.6p holds, and with f{X) = 
exp[— zS'(X)] for all z G C, we have the exact fluctuation symmetry 

(e-^^W) = (e-(i"-)S(^)) (3.11) 

with expectations in the nonequilibirum process starting from the equilibrium 
distribution po- 

Another way to get an exact fluctuation symmetry is to look back at p.9p 
with probabilities v = p = p equal to the stationary distribution of the nonequi- 
librium process. We then have from p.7p when combined with (|3.9p that in the 
nonequilibrium steady regime, for all T, 

{f{X)) = (e-'^(-^)-''(^°)+''(^^)/(eX)) (3.12) 

for irreversible entropy flux <j{X) = /? F{xt- , xt) and state function h{x) := 
logp(x) -I- pU{x). The exact symmetry (|3.12p invites to give special physical 
meaning also to that function h, but there is no convincing thermodynamic or 
operational meaning yet. That is also why asymptotic (in T f +oo) fluctuation 
symmetries have been more appreciated, obtained from p.l2p for / any positive 
function of (j(X), by taking the logarithm on both sides and using the bound- 
edness of the function h which makes it disappear when flnally dividing by T 
and letting T t +oo. 
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3.2 Relations with other aspects of nonequihbria 

The above techniques and relations are not new. Looking backward, it ap- 
pears that their main input has been the relation p.Sp . That has analogues 
for diffusion process [H [221 HI] , for dynamical systems [U [T71 dHl US] and also 
for non-Markovian processes [31 [301 [3T] as long as there is sufficient space-time 
locality to ensure a large deviation principle 0. The main origin of the fluc- 
tuation symmetry is therefore the identification of the entropy flux as marker 
of time-reversal breaking, [T71 [31 [51 [S] • We note next some other relations with 
aspects of nonequilibrium statistical mechanics. 

Quite some features of the close-to-equilibrium regime are easily deduced 
from the fluctuation symmetry. There are for example the Green-Kubo re- 
lations, with Onsager reciprocity as first explained in |32] following from an 
extended fiuctuation symmetry. The fluctuation-dissipation theorem with the 
Kubo formula [33| is a more general consequence. More globally, the validity of 
the McLennan ensemble close-to-equilibrium is another implication, see [331 ISi] • 

Wc will illustrate just one aspect which we have not seen stated as such, and 
which is useful. Start again from (|2.6p and take a function f{X) — g{QX)—g{X) 
in terms of another function g of interest. Then, 

{g{X)) = {g{GX)) + {{g{ex) g{X)) e'^W) (3.13) 

Imagine now that the action A in (|2.ip is small, so that the law P is just a small 
perturbation of the reference law and so that S = AQ — A is small. We can 
then expand the last term in p.l3p to bring 

{g{X)) = {gieX)) + {gieX) g{X)) {{g{QX) - g{X))S{X)r 

= {g{QX)) + {g{X)S{X)Y (3.14) 

where the last expectation, with the superscript (•)°, is with respect to the 
reference Pq and we have used that Pq is O— invariant. That linear order relation 
can be applied to the context of dynamical ensembles as we had it above, with 
@ time-reversal on trajectories X = (a;t,i S [0,T]). Take for example g{X) = 
0{xt) so that g{QX) = 0(a;o) for a state function O. We then obtain from 
p.l4p the Kubo formula 

{0{xt)) = (0(xo)) + {0{xt) S{X))l^ (3.15) 

where the expectations refer to the process P started from equilibrium po at 
time zero. Indeed, we should substitute in (|3.15p the expression p.Sp for S{X) 
with F = 0, /3t = P^at — a ~ £tO{i) and p, = v = po being the equilibrium 
distribution with potential U{x,a). Then, still using the flrst law (|3.2p . we 
arrive at the more familiar linear response expression 

(0(xt)) - (O(xo))^/ - (3.16) 
{0{xt) Six)):" = ^ dses^JO{x,)^U{x^,a)ro'' 
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Yet, it takes the combination p.8p - p.l5p to immediately understand why this 
formula is truthfully called fluctuation- dissipaiion relation. 

Moving beyond the linear response around equilibrium, makes it more dif- 
ficult to find specific consequences. Of course, the fluctuation relations hold 
unperturbed but there is no direct way to derive more specific results. In fact, 
it appears that one really needs more information about the time-symmetric 
part, T in (|2.4p . to move further, [HI [35], and that is also part of the motiva- 
tion of the next sections. 



4 Symmetry in dynamical activity 

We come to give examples of the fiuctuation symmetry (|2.7p . referred to in the 
title of the paper as active because they deal with the dynamical activity. 

4.1 Boundary driven Kawasaki dynamics 

Take K ~ {0, where states are particle configurations x — {x{i), i G 

{1, 2 . . . , £}), x{i) = 0, 1, interpreted as vacant versus occupied sites on a lattice 
interval. In other words we are speaking about indistinguishable particles sub- 
ject to exclusion on a lattice interval. The dynamics has two parts. First, bulk 
exchange of neighboring occupations, for inverse temperature /? > 0, 

k{x, y) = cxp - V{x)] 

when y{j) = x{j) for all j except for y{i) ~ x(i + + 1) = x(i) for some 

i = 1,2. — 1. The interaction between neighboring sites is ruled by the 
potential 



U{x) = X! + 1) 



where k S M is the coupling parameter. Apart from that interacting diffusion 
part to the dynamics, there arc also the reactions at the boundary sites; there 
is creation and annihilation of particles at i = 1 and i = L, with rates 

k{x, y) = exp - ^ [U [y) - U{x)] exp ^ y) 

for y{j) — x{j) except for j = 1 where y{\) — 1 — a;(l), or for y[j) — x{j) except 
for j ~ L where y{L) = 1 — x(L), and 

A{x,y) = +{a + d)wheny{l)^l,x{l)=0,y{j)^x{j),j:f^l, 

= -{a + 6) when 2/(1) = 0, x{l) = 1, yO') = a;(j), j ^ 1, 

= +(a - 5) when y{L) = \,x{L) = 0, y{i) = x{j),j ^ L, 

= -{a -6) when y{L)^0,x{L)^ I, y{j) ^ xij),j (4.1) 
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for some fixed parameters a, 5 £ R. The interpretation is that there are left and 
right chemical potentials /if,/ir of the particle reservoirs of the interval with 
fiifi :— a + 5, firp :— a — 6. For all other transitions k{x, y) = 0. As a result, 

k{x,y) = koix,y) exp[- J{x,y)] 

with ko{x,y) = exp[S{y)^S{x)]/2,S{x) :^ -l3U{x)+aN{x), N{x) := Eti^i^) 
(number of particles in the system for state x), and current 

J{x,y) = +1 when a particle enters at ?' = 1 

= —1 when a particle leaves at i = 1 

= +1 when a particle leaves at i = L 

= —1 when a particle enters a.t i ~ L (4-2) 

and zero otherwise. In other words, J{x,y) = Jr{x,y) — Ji{x,y) with Ji{x,y) 
the current of particles into the left reservoir under the transition x ^ y, and 
Jr{x,y) the current of particles into the right reservoir. 



For (5 = (and only for 5 = 0) there is detailed balance with grand-canonical 
ensemble ^ 

Po{x) = — cxp5(a;). 

Then, a/ (3 is the chemical potential of both particle reservoirs left and right. 
That equilibrium process determines our reference distribution Pq. Noncqui- 
librium arises from taking 5 7^ 0, which makes the chemical potentials in the 
imagined left and right particle reservoirs different. We can start the noncqui- 
librium process from the same po, giving our distribution P, but asymptotically 
in time a noncquilibrium steady regime will develop. In particular it can be 
proven that for 5 > there will be a steady particle current from left to right. 
See for example |37] for the details of the standard fluctuation symmetry as in 
the previous section. 



The decomposition (|2.4I) here gives 

S{X) = 5[J,{X)- MX)] (4.3) 

with S{QX) = —S{X) for 6 time-reversal, and Ji{X) := Jt{xt- , Xt) the net 
number of particles that escape from the interval to the left particle reservoir. 
Note that Jr[X) + Je{X) = —M{xt) + A/'(xo), the change of the number of 
particles in the system. 

For the time-symmetric part of the action we have computed from p.lOp that 

r{X) = 2j At[Bi{xt;a,5) + BL{xua,5)] (4.4) 
Jo 
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where 

+ (e(-+«)/2_e'^/2)(g/^K/2_^)^(2) 

+ {(e-'3«/2_l)(e-(a+^)/2_g-a/2) 

- (e^''/^ - l)(e('^+'')/2 _ e'^/2)}x(l)a;(2) 

and 

Bl{x; a, S) := e^^-*)/^ - e"/^ + {e-(--«)/2 _ e(-*)/2 + ea/2 _ e-a/2|^(^) 

- e'^/2)}a;(L - l)x(L) 

We next apply the mirror symmetry F through which {TX)t{i) = Xt{L — i + 1). 
Observe that in that mirror Jf,{X) — Jr{TX), ST ~ —S. We can thus compute 

R{x) = \ (r(rx) - r{x)) = At r{xt) (4.5) 

from the expected difference in transitions (jumps in and out of the system) left 
versus right, to find 

r(x) = -2 sinh ^ ((e^'^'^/^ - l)e-''^^ + (e'^'^/^ - l)e"/2^ {x{L)x{L - 1) - x(l)a;(2)) 
+ 2( sinh -sinh ^^-^)(a;(L) -a;(l)) 

+ 2e''/^ sinh^ie^'^/^ -l){x{L-l)-xi2)) (4.6) 

which is of course also odd in the driving field 5. For the boundary driven 
symmetric exclusion process we must take the coupling k = 0, and only survives 
the term 

r«=° = 2( sinh - sinh ^i^)(x(L) - .t(1)), (4.7) 

given entirely in terms of the difference in occupations at the outer sites. 

It follows from the general analysis in Section [2] that R{X) in (|4.5p verifies 
the fluctuation symmetries p.7p - p.l0p . which is a non-trivial general identity 
whose meaning refers to the reflection-antisymmetric part in the dynamical 
activity (|4.4p . In particular, that identity (|2.7p for that same R in (|4.5P " (|4.6p 
remains strictly valid even when modifying the interaction potential U in the 
bulk of the system. On the other hand, applying the general consequence that 
{R{X)) > 0, or J2x ''(^') Pi^) — 0: to (|4.7p only gives the well known fact that 
the density is larger at the side of the largest chemical potential. 



12 



4.2 Boundary driven zero range process 



We discuss next the application of fluctuation symmetries to a bosonic version 
of the previous example, where particles diffuse without exclusion principle. 
Consider again a one-dimensional channel composed of L cells in which we 
observe occupation numbers n{k) e N, k ~ 1,. . . ,L. The particle configuration 
X = (n(l), . . . , n(L)) can change in two ways. First, via bulk hopping x — >■ 
x~ei + ei±i at a rate w{n{i)), where stands for the particle configuration with 
one particle in cell i and zero elsewhere. The choice w{n{i)) oc n(i) corresponds 
to independent particles. Secondly, at the boundaries, the channel is connected 
to particle reservoirs with chemical potentials /xi = /i, /x^ = /i, respectively. The 
transition rates for the creation/annihilation of particles at the boundary sites 
arc then 

k{x,x — ei) = siw{ni) 
k{x,x + ei) = f 1 := ri e'' 
k{x,x — eL) — SLw{nL) 

k{x,x + eL) = fL-^TLe-^ (4.8) 

The rates for these transitions evoke the chemical potentials at the boundary 
walls from = \og{fL/ sl) and jl = log(fi/si). We assume that si/ri = sl/tl 
so that, for (5 = 0, we have the equilibrium situation where the chemical po- 
tentials left and right become equal. Of course wc could have chosen also to 
modify the si,sl, but it appears physically most accessible to change the in- 
coming rates ri — ?> fi,ri ^ to achieve a noncquilibrium regime as wc also 
do in the next section. In fact, to make the equilibrium left/right symmetric we 
also take si = SL,ri = r^. 

The corresponding stationary distributions po (at ^ = 0) and p (at general 5) 
arc product distributions but that will not be used in the following. 



We consider trajectories X = {xt,t G [0,T]). Both the equilibrium Po and 
the noncquilibrium process P start from the same equilibrium distribution po- 
The action (|2.ip is easily calculated to be 

A{X) = log ^ li{X) + log ^ Il{X) -T {n+TL- n - rl) (4.9) 

where e.g. /f,(^) indicates the number of particles entering the system from 
the left reservoir for the path X. As wc apply time-reversal 0, wc obtain the 
time anti-symmetric part of the action S{X) = A{QX) — A{X) 

S = log- (It- It) + log— {r ~ II) 

^6{r-Je) (4.10) 

where e.g. Je :— li_ — is the net number of particles that have escaped 
to the left particle reservoir during [0,T]. As usual and as explained before. 
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that entropy production satisfies the exact fluctuation symmetry (|2.9p . For the 
asymptotic form, one must be more careful because of the unbounded number 
of particles; see [55] • Here we are however more interested in the dynamical 
activity. 

That is given in the time-symmetric term T(X) = A{QX) + A{X), 

T=-d {It +lt)+d {II +11) -2 {h +fL-rL- ri) T (4.11) 

which is the analogue to (j4.4p . As there, we now apply the mirror transformation 
r, reversing left /right. First note that again S is antisymmetric under F, ST = 
—S. On the other hand, we have 

T{TX) - T{X) = 25 {it +li-r_,- II) (4.12) 

exactly proportional to the difference in dynamical activity 

A(X) := Z+Il-li- It 

between the right and left boundary. Following the logic of (|2.7I) . that suffices 
for TT — T oc A to satisfy a fluctuation symmetry (|2.10p up to a total time- 
difference: when fQ = / is time-symmetric, then 

(/(FX)) = (/(X)e*^m) (4.13) 

for all times T, where we start the nonequilibrium process at time zero from po- 
For example taking / = A, to first order in S, 

{A{X)} = -^-{A^{X)}l'' (4.14) 

which is formally similar to a Green-Kubo relation but the observable A is time- 
symmetric. 

It is in fact true for all (5 > that (A) < which means that the greatest 
activity is to be found at the boundary side of the largest chemical potential. In 
other words, as for the boundary driven Kawasaki dynamics also for zero range, 
the particle current can be said to be directed away from the region of largest 
activity. These statements all hold for any form of the bulk rate w and are 
quite independent of the usual statements involving the fiuctuation symmetry 
of entropy production or of currents. 

5 Spinning Lorentz Gas 

The Spinning Lorentz gas (SLG) is a classical mechanical model of particle 
scattering in 2D; it is actually an interacting version of the normal Lorentz gas 
[39] , which is a well known example of deterministic particle diffusion [40] [41] . 
The SLG has the additional feature of providing local thermalization of the 
wandering particles along with the scatterers; a complete description of this 
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Figure 1: In the Spinning Lorentz Gas (SLG) M disks with radius one and center 
fixed in a triangular lattice rotate freely and exchange energy with particles (of 
mass one) via elastic collisions [42] . The particles move mechanically inside the 
slab of length L with periodic boundary conditions in the vertical coordinate. 
The slab is placed among thermo-chemical reservoirs (ideal gases) with (for the 
present paper) equal inverse temperatures f3 and different chemical potentials 
^J'^=l,L■ The particles can enter and leave to/from the reservoirs at the left and 
right boundaries. 



and the coupled energy and mass transport properties of the SLG model can 
be found in [32] • As a matter of fact, the validity of the fluctuation theorem 
for the entropy production and for the joint distribution of currents has been 
tested for this model before, of course taking into account the limitations due 
to the unbounded kinetic energy, see [43j . Also, a precise meaning to the state 
function h of p.l2|) can be found in the SLG context, for the exact fluctuation 
symmetry case |43] . Specifically, by including such model here we aim to aid 
extending the study of the nonequilibrium fiuctuation relations within the realm 
of the time-symmetric variables (the dynamical activity). For this purpose we 
use closely the relations obtained for the previous example in Section 14.21 to 
apply them in this somewhat more realistic transport model. 



As illustrated in Figure [U in the SLG the array of scatterers is connected 
to thermo-chemical reservoirs, with nominal chemical potential = 1,L and 
with inverse temperatures /3. This setting drives the system into a nonequi- 
librium stationary regime, when the chemical potentials of the reservoirs are 
different. 

The SLG model is of course more mechanical than the boundary driven zero 
range model of the previous section. We want to connect them however. At 
the walls, the particles can enter and can leave the system. The rates at which 
particles enter are in general related to the mean density u of their reservoir as 
oc see also [44]. In the nonequilibrium setting, there is a reservoir chemical 

potential difference, given by /3A// = /3(a*l — /"i) = log(itL/wi); and hence, in 
the notation of the previous section, wc put = ^^i e'^'^^, or 2S = —(SAji. 
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Spinning Lorentz Gas (NEQ). L ^ 40; l3An = 0.20 



3 
2.5 



Spii.i.iug Lurtiili! r.aa (NEQ). L = 40^ /3A(i = 0.20 
is-' = .10 ' 
/3-' = Kill 



= 50 
identity 
3"^ = 100 




Figure 2: The fluctuation symmetry for the dynamical activity is tested numer- 
ically in nonequihbrium simulations of the SLG. In the inset, the probability 
distribution Pt (R) of R in equation (j5.ip is given. In the vertical axis we have 
the calculation of the functional 15.21 The slab length is L = 40, with reservoir 
chemical potential difference /3A/i = 0.2, and reservoir temperatures /3~^ = 50 
(crosses) and {3^^ = 100 (stars) giving identical result. 



We now wish to conjecture that identical fluctuation relations as (|2.7p - (|2.10p 
hold for the dynamical activity as we had it for the boundary driven zero range 
process before, in particular in the version (|4.13p . One therefore looks back at 
the expression (|4.12[) . More precisely, we take the time-symmetric variable 



R 



2T 



{{It + It) - {II + O) 



(5.1) 



Figures [5] and [3] show the validity of the time-symmetric fluctuation theorem 
for the probability Pt {R) in the SLG model. This distribution is obtained via 
molecular dynamics simulations of the system in a nonequihbrium stationary 
state, with different reservoir chemical potentials. In these figures we plot the 
functional 

1 . Pt{R) 



Ut {R) 



In- , 

T Pt{- 



R) 



(5.2) 



The measuring time was a large value of T = 4.0, which means that in 15.21 one 
understands the flucutation symmetry in the asyptotic sense; in the same time 
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units, the average time between collisions in the gas is ^ 2.5 x 10"'^. In the 
first case (Fig. [2]) stationary nonequihbrium is obtained by a chemical potential 
difference /SAfi = 0.20 and for two different temperatures. 

The second case (Fig. ^ corresponds to a larger driving f3Afi ~ —0.45; this 
gives a fluctuation theorem interval in which the distribution is visibly non- 
Gaussian. 

As in the remark around (|4.7p , the dynamical activity in (jS.ip is proportional 
to the number of transitions at the wall; in other words, it is proportional to 
the local boundary density. Since the temperature in this case is uniform, the 
activity fluctuations are simply related to density fluctuations in the stationary 
profiles. Thus, when measuring the differences in dynamical activity in (|5.ip 
one obtains asymmetric statistics due to the asymmetry in the nonequihbrium 
density profile along the slab in fig. (TJ by the condition set in the reservoirs. 
Thus, one observes that the variable in (j5.1|) is sensitive to the breaking of spatial 
symmetry, which is another valid aspect of nonequihbrium; this is related but 
also complementary to the breaking of time-reversal symmetry. That is then 
the physical meaning of the involution F introduced in Section [2j as previously 
noted in [T5] . 



6 Conclusions 

Fluctuation relationos in nonequihbrium extend to the realm of time-symmetric 
observables. These observables are however antisymmetric for other transfor- 
mations such as spatial symmetries, the breaking of which is responsible for 
the time-reversal breaking. The inhomogeneity in dynamical activity exactly 
picks up that (other) symmetry breaking. In particular, we have checked via 
computer simulations that there is a fluctuation symmetry for the difference 
in dynamical activity at the boundaries of the spinning Lorcntz gas in station- 
ary nonequihbrium. These were guessed from a mathematical analysis of the 
boundary driven Kawasaki and zero range process. We would like to call these 
active fluctuation symmetries as they do not involve the dissipative currents, but 
rather that complementary aspect of frenesy or dynamical activity. As we have 
seen, the logic and mathematical derivation of these additional time-symmetric 
fluctuation symmetries is completely analogous to the existing fluctuation re- 
lations for dissipated work and entropy fluxes, except that an additional sym- 
metry is involved; most simply a mirror or reflection symmetry, which basically 
is equivalent to reversing the driving field. The very fact that reversing time 
is not equivalent with reversing the driving field is responsible for non-trivial 
nonequihbrium behavior beyond the linear regime. 
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Figure 3: The validation of the time-symmetric fluctuation theorem for the 
case of non-Gaussian fluctuations of the dynamical activity difference R, in the 
SLG in stationary nonequilibrium. The chemical potential difference in the 
reservoirs is /3A/x = —0.45, with /3 = 1/150 and slab length L = 40. The inset 
shows the probabihty Pt{R), measured from the computer simulation with a 
large measuring time T = 4.0. The interval of fluctuations around zero is far 
from the average value, where one distinguishes non-Gaussian behavior. In the 
main plot, the crosses show the evaluation of the fluctuation theorem for the 
probabilities in the inset; these data fit to a straight line with slope close to one, 
m = 0.99735 ±0.01265. 



18 



References 



[1] D.J. Evans, E. G. D. Cohen and G.P. Morris, Probability of second law 
violations in steady flows. Phys. Rev. Lett. 71, 2401-2404 (1993). 

[2] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995); J. 
Sat. Phys. 80, 931 (1995). 

[3] C. Maes, The fluctuation theorem as a Gibbs property. J. Stat. Phys. 95, 
367-392 (1999). 

[4] J.L Lebowitz and H. Spohn, The Gallavotti-Cohen fluctuation theorem for 
stochastic dynamics. J. Stat. Phys. 95, 333 (1999). 

[5] G. E. Crooks, Entropy production fluctuation theorem and the nonequilib- 
rium work relation for free energy differences. Phys. Rev. E 60, 2721-2726 
(1999). 

[6] C. Maes and K. Netocny, Time-reversal and Entropy. J. Stat. Phys. 110, 
269-310 (2003). 

[7] C. Maes, On the origin and the use of fluctuation relations for the en- 
tropy. Seminaire Poincare 2, 29-62 Eds. J. Dalibard, B. Duplantier and V. 
Rivasseau, Birkauser (Basel), 2003. 

[8] P.G. Bergman and J.L. Lebowitz, New Approach to Nonequilibrium Pro- 
cess. Phys. Rev. 99, 578-587, 1955. 

[9] S. Katz, J.L. Lebowitz, and H. Spohn, Stationary nonequilibrium states for 
stochastic lattice gas models of ionic superconductors. J. Stat. Phys. 34, 
497-538 (1984). 

[10] T. Harada and S.-Y. Sasa: Equality connecting energy dissipation with vi- 
olation of fluctuation-response relation, Phys. Rev. Lett. 95, 130602 (2005). 

[11] B. Derrida, Non-equilibrium steady states: fluctuations and large devia- 
tions of the density and of the current. J. Stat. Mech. P07023 (2007). 

[12] H. Tasaki, Two theorems that relate discrete stochastic processes to micro- 
scopic mechanics. larXiv:0706.1032k ^l [cond-mat.stat-mech]. 

[13] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997); Phys. Rev. E 56, 5018 
(1997). 

[14] M. Baiesi, E. Boksenbojm, C. Maes and B. Wynants, Nonequilibrium lin- 
ear response for Markov dynamics, I: jump processes and overdamped dif- 
fusions, J. Stat .Phys. 137, 1094-1116 (2009). 

[15] J. P. Garrahan, R.L. Jack, V. Lecomte, E. Pitard, K. van Duijvendijk, 
and F. van Wijland: First-order dynamical phase transition in models of 
glasses: an approach based on ensembles of histories, J. Phys. A: Math. 
Gen. 42, 075007 (2009). 



19 



[16] R. Jack, J. P. Garrahan, D. Chandler: Space-time thermodynamics and 
subsystem observables in kinetically constrained models of glassy materials, 
J. Chem. Phys. 125, 184509 (2006). 

[17] F. Bonetto, G. Gallavotti, and P. Garrido, Chaotic principle: An experi- 
mental test. Physica D 105, 226 (1997). 

[18] C. Macs and M.H. van Wicrcn, Time-symmetric fluctuations in nonequi- 
librium systems. Phys. Rev. Lett. 96, 240601 (2006). 

[19] W. De Roeck and C. Maes, Symmetries of the ratchet current. Phys. Rev. 
E 76, 051117 (2007). 

[20] C. Maes, K. Netocny and M. Verschuere, Heat Conduction Networks. J. 
Stat. Phys. 111,1219-1244 (2003). 

[21] J.W. Silverstein and P.L. Combettes, Signal detection via spectral theory 
of large dimensional random matrices. IEEE Transactions on Signal Pro- 
cessing 40(8), 2100-2105 (1992). 

[22] M. Colangeli, C. Maes and B. Wynants, A meaningful expansion around 
detailed balance. J. Phys. A: Math. Theor. 44, 095001 (2011). 

[23] One should have in mind that Q is time-reversal on trajectories and F is 
a spatial transformation. In many cases, the breaking of time-symmetry is 
directly related to the breaking of a spatial symmetry. 

[24] K. Sekimoto, Stochastic Energetics. Lecture Notes in Physics 799, Springer 
(2010). 

[25] F. Ritort, C. Bustamentc and I. Tinoco,Jr., A two-state kinetic model for 
the unfolding of single molecules by mechanical force. Proc. Nat. Ac. Scienc. 
99, 13544-13538 (2002). 

[26] J. Kurchan, Fluctuation theorem for stochastic dynamics. J. Phys. A: Math. 
Gen. 31, 3719-3729 (1998). 

[27] C. Maes, F. Redig and A. Van Moffaert, On the definition of entropy pro- 
duction, via examples. J. Math. Phys. 41, 1528-1554 (2000). 

[28] D. Ruelle, Smooth Dynamics and New Theoretical Ideas in Nonequilibrium 
Statistical Mechanics. J. Stat. Phys. 95, 393-468 (1999). 

[29] C. Maes and E. Verbitskyi, Large deviations and a fluctuation symmetry 
for chaotic homeomorphisms. Comm. Math. Phys. 233, 137-151 (2003). 

[30] D. Andricux and P. Gaspard, The fluctuation theorem for currents in semi- 
Markov processes. J. Stat. Mech.:Theory and Experiment, P11007 (2008). 

[31] C. Maes, K. Netocny and B. Wynants, Dynamical fluctuations for semi- 
Markov processes. J. Phys. A: Math. Theor. 42, 365002 (2009). 



20 



[32] G. Gallavotti, Extension of Onsagers Reciprocity to Large Fields and the 
Chaotic Hypothesis. Phys. Rev. Lett. 77, 4334-4337 (1996). 

[33] R. Kubo, The fluctuation-dissipation theorem. Rep. Prog. Phys. 29, 255- 
284 (1966). 

[34] T. S. Komatsu and N. Nakagawa, An expression for stationary distribution 
in nonequilibrium steady states. Phys. Rev. Lett. 100, 030601 (2008). 

[35] C. Maes and K. Netocny, Rigorous meaning of McLennan ensembles. J. 
Math. Phys. 51, 015219 (2010). 

[36] C. Macs, K. Netocny, and B. Wynants, On and beyond entropy production; 
the case of Markov jump processes. Markov Proc. Rel. Fields. 14, 445-464 
(2008). 

[37] C. Maes, K. Netocny and B. Shergelashvili, A selection of nonequilibrium 
issues. In: Methods of Contemporary Mathematical Statistical Physics, 
Ed. Roman Kotecky, Lecture Notes in Mathematics 1970, pp. 247-306, 
Springer, 2009. 

[38] R. J. Harris, A. Rakos and G.M.Schiitz, Breakdown of Gallavotti-Cohen 
symmetry for stochastic dynamics. Europhys. Lett. 75, 227 (2006). 

[39] H. A. Lorentz, Proc. Amst. Acad. 438 (1905). 

[40] P. Gaspard, Chaos, Scattering, and Statistical Mechanics (Cambridge Uni- 
versity Press, Cambridge, 1998). 

[41] Hai^d Ball Systems and the Lorentz Gas, edited by D. Szasz (Springer Ver- 
lag, Berlin, 2000). 

[42] H. Larralde, F. Leyvraz and C. MejaMonasterio, Transport properties of a 
modified Lorentz gas. J. Stat. Phys. 113, 197-231 (2003). 

[43] A. Salazar, F. Leyvraz and H. Larralde, Fluctuation theorem for currents 
in the Spinning Lorentz Gas. Physica A 388, 4679-4694 (2009). 

[44] F. Reif, Fundamentals of Statistical and Thermal Physics (International 
Student Edition), McGraw-Hill, Tokyo, 1965; pp. 272. 



21 



